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I. INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

The XV-15 is an advanced tilt rotor research aircraft (Figure 
1.1) currently under development by the National Aeronautics ana 
Space Administration and the U.S. Army Air Mobility R§D Labora- 
tory. An important element of the aircraft development is an ex- 
tensive full-scale wind tunnel test to be conducted in the Ames 
Research Center 40- by 80- foot wind tunnel. Since previous wind 
tunnel scale model and full-scale tests of the tilt rotor concept 
have been effective in previous years, it is desired to conduct the 
forthcoming XV-15 tests in a comprehensive manner to continue to 
minimize uncertainties in system characteristics. 



Figure 1.1 The NASA/ Army XV-15 
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This study is a determination of specific test requirements 
which impact the conduct of these tests. This research is based on 
a mathematical model of the XV-15. This model is used to analyze 
the dynamic characteristics of the vehicle at 190 knots in the 
wind tunnel, 

1.2 METHOD OF APPROACH 

The method of approach used for this work is schematically 
shown in Figure 1.2. The mathematical model was converted into a 
simulation and a state vector format. These two reformulations 
were then combined with existing data on the input control channels 
and instrumentation characteristics to provide a model for the test 
system. 

The next steps were the evaluation of the test model to deter- 
mine basic requirements on the data analysis algorithms. In par- 
ticular, this step produced an evaluation of mode identif iability 
which set further requirements on instrumentation and inputs. 

Originally, the effort was directed to a basic rotor/cantilever 
wing model. During the course of this effort, coordination with 
the test agency indicated the need for results which expanded this 
method of approach to coupled modes between both wings. The steps 
were repeated for this case. 

1.3 PRINCIPAL CONTRIBUTIONS AND CONCLUSIONS 

The following analytical tools have been developed: 

(1) An XV-15 digital simulation incorpoi ating a simplified 
tunnel support model, control system loop, measurement 
lags, gust disturbances, and sensor noise. Time histor- 
ies are generated from this simulation which provide a 
correlation base for the tunnel tests, as well as a means 
of evaluating various data processing methods. 
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Figure 1.2 Study of XV-IS Tunnel Test Procedures and Analysis 




(2) Specialization of existing data analysis programs to the 
high order XV- IS dynamical model. The programs so speci- 
alized at various stages of this work are: (a) the trans- 
fer function program, (b) a time series analysis program 
(FFT, auto-correlation, transfer functions), and (c) SCIDNT, 
an advanced maximum likelihood parameter identification 
program. 

(3) Several auxiliary programs have been developed to provide 
estimates of damping from transfer functions, as well as 
calculation of modal decomposition of system response 

(to identify modes), particularly useful for analysis of 
coupled wing responses. 

Application of these programs to the basic symmetric, nine 
degree-of- freedom mathematical model have produced the following 
conclusions (see Table 1.1): 

(1) If the test configuration at 190 knots is such that con- 
trol inputs are limited to collective pitch and flaperon 
at frequencies below approximately 5 Hz and measurements 
are made only of the wing bending accelerations, and 
q^ 2 » then sufficient information for calculating modal 
frequency and damping accurately is available only for 
the q^j_ and q ^2 >nodes. (Necessary input amplitudes to 
achieve a desired noise-to-signal ratio for different 
levels of wind gust and speed of actuator response can 
be obtained from Figures 3. 3-3. 7 for the cantilever wing 
and Figures 4.6-4. 8 for the coupled wing.) 

(2) If the above test configuration could be improved to 
allow flaperon input frequencies in the vicinity of 10 Hz, 
then information would be sufficient to compute modal 
frequency and damping of the wing torsion mode from 
measurements of the wing bending accelerations. 
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Table 1.1 

Frequency and Damping Information Available in Various Test Configurations for 

the Nine Degree-of- Freedom Symmetric Model 
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(3) Similarly, collective pitch input frequencies near 20 Hz 
would provide information for the frequency and damping 
of the rotor coning mode, 6. 

(4) If instrumentation is provided to measure the wing tor- 

sion degree-of- freedom, p, and if cyclic pitch control at 
18.6 Hz is possible, then the frequency and damping of 
the upper rotor inplane bending, could be determined. 

(5) Even under the best of circumstances (i.e., any of the 

four controls could be excited at any frequency and 
measurements could be taken of all nine degrees-of -freedom) 
sufficient information is not present in the transfer 
functions to calculate the modal frequencies and dampings 
of the upper and lower gimbal modes, and B.j^, and the 
lower rotor inplane bending mode, It is because 

these modes are well damped that there are no resonant 
peaks associated with them in the transfer functions. 

Consideration of the coupled wing motions in the presence of 
control and measurement calibration errors produced the following 
conclusions : 

(1) The test configuration of exciting one wing and measur- 
ing the response of that wing makes the determination of 
frequency and damping of individual modes extremely 
difficult since the symmetric anv. asymmetric modes 
usually result in a single resonant peak in the frequency 
response. 

(2) Symmetric and asymmetric motions can be separated if 
either both wings are excited and/or the responses of 
both wings are measured. 


6 



*»■ 


(3) Errors in calculating a mode's frequency and damping 
from its resonant peak may still occur if the motion of 
other modes is significant. 

(4) Measurement and control calibration errors less than 20% 
do not produce significant errors in calculating damping 
values compared to the problem stated in (3) above. 

1.4 SUMMARY OF REPORT 

The objective of this report is to present principal modeling 
and analysis approaches and reference material for tunnel test 
planning. Chapter II discusses the modeling of the aircraft, 
instrumentation, and controls. Chapter III reviews results of the 
rotor/cantilever wing model, and Chapter IV presents the coupled 
wing results. Chapter V gives some examples of data prediction 
with system identification techniques. Chapter VI presents de- 
tailed conclusions and recommendations. 

The appendices provide the detailed analysis and figures upon 
which the text draws. Appendix A is the derivation of the support 
equations. Appendix B presents the principal derivations for 
transfer function analysis. Appendix C is the plots of transfer 
functions . 
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II. TILTING PROPROTOR MATHEMATICAL MODEL 




An analytical evaluation of tilting proprotor test require* 
ments is based on a comprehensive description of the dynamics and 
aerodynamics of the subject vehicle. This section discusses the 
basis of that description (the simulation, detailed in [2]j, and 
modifications to this simulation which approximate typical test con- . 
figuration constraints and disturbances. Section 2.1 reviews essen- 
tial elements of the basic rotor/wing model. Section 2.2 discusses 
the influence of the aircraft tunnel support restraint. Coupled wing 
interactions are approximated in Section 2.3, followed by Section 
2.4, outlining the gust simulation method. Section 2.5, modeling of 
instrumentation, and Section 2.6, control system representations. 

2.1 REVIEW OF BASIC ROTOR WI.NG MATHEMATICAL .MODEL 

The basic mathematical model of the XV-IS aircraft dynamics 
has been established by Johnson [2-5]. A conclusion of Reference 5 
was that, in general, the basic rotor dynamics were satisfactorily 
described by a nine degree-of-freedora model (the first bending 
mode and rigid pitch mode of each blade, gimbal pitch and roll 
angles, and rotor speed perturbation), but that in some cases it 
may be reduced to six degree-of -freedom by using the quasi-static 
blade-torsion approximation, which is discussed in Reference 2. 

For the purposes of this study, it was determined that the six 
degree-of-freedom model was satisfactory. Hence, the complete rotor/ 
pylon/wing model consists of the six rotor and pylon degrees-of- 
freedom described above, plus three degrees -of- freedom for the wing, 
vhich are lowest mode vertical bending, chordwise bending, and 
torsion. These nine degrees -of-freedom are listed in Table 2.1 for 
jase of reference. 


PAGE''lNTEN.i.,:a_y _ ; :;c 
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Table 2.1 

Rotor, Pylon, Wing Degrees-of-Freedom and Controls 


SYMBOL 

DESCRIPTION 


Degrees -of - Freedom 

01c 

Longitudinal component In nonrotating frame of lowest blade 
bending mode 

»ls 

Lateral component in nonrotating frame of lowest blade bending 
mode 


Coning of the rotor 

®GC 

Gimbal pitch angle 


Gimbal roll angle 


Rotor rotational azimutL perturbation 


Wing vertical bending 

^2 

Wing chordwise bending 

P 

Wing torsion 


Controls 

®0 

Blade collective pitch 

9s 

Blade longitudinal cyclic pitch 

9c 

Blade lateral cyclic pitch 


Wing flaperon deflection 



Several aspects of this model deserve further comment. The 
lowest blade bending modes are essentially inplane bending (lead-lag) 
during high speed, axial flow. The effect of collective lag of the 
blades appears in the rotor speed perturbation degrees -of -freedom. 

The effects of engine and transmission dynamics are included in the 
equations of motion, but they do not introduce additional degrees - 
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o£-£reedoia. The equations o£ motion include the e££ects o£ gusts 
and £our control variables (see Table 2.1): collective, longitudin- 
al cyclic, and lateral cyclic pitch o£ the blades and wing flaperon 
deflection. 

The computer program developed by Johnson [3] calculates the 
equations of motion in the general form 

^2^ + A^x + A^x « Bu + B^g (2.1) 


where 

X is the state vector composed of either the nine sym- 
metric or asymmetric degrees-of-freedom, 

u is the vector of control variables, 

g is the gust vector, 

A2 »Aj^,Ajj are the ’’mass”, ’’damper”, and ’’spring” matrices, 
respectively, and 

B,Bq distribute the controls and gusts, respectively, 
among the states. 


This equation is readily converted to the standard state-space 
form 


X ■ Fx + Gu + Tg 


where 


X 



F 




I 


-A 


-1 

2 


A 


1 


( 2 . 2 ) 


(2.3) 


(2.4) 
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Measurements, y, of the degrees-of-freedom are modeled as lin- 
ear combinations of the state and control: 

y ■ Hx + Du + D"‘g (2.7) 

Except where stated otherwise, it is assumed that the accelerations 
of the wing bending degrees-of-freedom, and are measured by 
accelerometers and that direct measurements of other degrees-of- 
freedom are available as needed for the sake of discussion. 

The frequency response of the measurements to the controls or 
gusts are available from the Laplace transfer functions 


^ - H(sl-F)'^G ♦ D 

(2.8) 

Xjlj - H(sl-F)‘^r • D' 

(2.9) 


which follow directly from Eq. (2.2). The results obtained 
from these functions are discussed in Chapters III and IV. 

2.2 EFFECT OF SUPPORT FLEXIBILITY 

A preliminary requirement for this study was to determine if 
the flexibility of the support structure in the 40 x 80 wind tunnel 
would introduce modal frequencies near frequencies of interest. If 
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so, the isolation of the modal frequencies of interest and the 
I determination of the modal damping factor would be made more diffi- 
cult. 


Displacement of the support attachment points in a horizontal 
plane was considered. Vertical displacement was not considered 
because the support is far more rigid in that direction. Thus, the 
basic equations of motion were expanded to include two more degrees - - 
of-freedom: longitudinal and lateral translation of the aircraft, 

Xg and respectively. Aircraft yaw is also possible, but data 
were not available to adequately specify the frequency and damping 
parameters for the yaw degree-of -freedom, so it was omitted from 
the study. 


The support equations will be only briefly described here; a 
detailed derivation of the equations appears in .Appendix .A. The 
direct aerodynamic effects on the x^ and y^ degree-of -freedom are 
negligible compared to the structural kinetic effects; therefore. 


only the latter 

are retained in 

the equations. 

The result may be 

written simply 

as 




A 

• 

A 



s s 



- F 

(2.10) 


s 

s 
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* 

^ • 

* 



m 9 
s ^s 

+ C y 

♦ K y 


(2.11) 


where the asterisk superscripts denote nondimensional quantities. 

The mass coefficient, m , is estimated from a weight analysis 

of the XV-IS. The terms F and F are composed of all the 

*s ^s 

cross-coupling terms between the support and the remaining degrees- 
of-freedom, plus the effects of controls and gusts, all of which 
were available from the existing computer program. Therefore, only 
the damping coefficients, C and C , and the spring coefficients. 


and K , 

' c 


S s 

remain to be specified. This was done by assuming 2% 
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structural damping and specifying the frequency of the support 
oscillations with the rotor off. These frequencies were taken as 
3.5 Hz in the longitudinal direction and 4.0 Hz in the lateral, 
based on earlier experience in the 40 x 80 tunnel [6]. 

The result is shown in Figure 2.1, which compares the pole 
locations of the symmetric nine degree-of-freedom model to the same 
model with the support degrees-of 'freedom added. As can be seen, 
the effect is an increase in the modal frequencies of the rotor- 
coning CSq) mode (from 2.6 to 2.9 per rev); of the upper blade lag 
mode (from 2.4 to 2.6 per rev); and of the wing chordwise 
bending Cq^2^ (from 0.67 to 0.84 per rev). The frequencies 

of the remaining modes are unchanged. The two new modes associated 
with Xg and y^ are at low frequencies (0.12 and 0.15 per rev, 
respectively) and are not expected to cause any difficulties in 
identifying the frequencies and damping of other modes. Also, even 
though the S^, and q^ 2 * frequencies are altered, there 

still remains adequate separation of frequencies to permit discrim- 
ination of individual modes. Therefore, the conclusion is that 
the support flexibility does not "change the basic characteristics 
of the XY-15 structural damping or frequencies and need not be 
considered in the remainder of this study. 

2.3 MODELING OF COUPLED WING RESPONSES 

A major part of this study was to ascertain the effects of 
exciting one or both rotors (or wings) , the effects of measuring 
the iresponse of one or both wings, and the effects of calibration 
errors in the measurements of the control inputs and the wing 
responses on the ability to correctly compute the frequency and 
damping of important modes. The nine degrees-of-freedom of the 
right rotor/pylon/wing system are coupled to the nine degree '-of- 
freedom of the left system, principally by means of the transmission 
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Figure 2.1 Comparison of Pole Locations for Symmetric Equations 
With and Without Support Degrees -of - Freedom 
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interconnect shaft. This left and right system of eighteen degrees- 
of-freedom may be decoupled into symmetric and asymmetric systems 
of nine degrees -of-freedom each. This procedure is advantageous 
because it aids the understanding of the complete system, simplifies 
the mathematics, and reduces the computational effort. 

There are nine cases to be investigated (Table 2.2). It is 
assumed that the symmetric and asymmetric responses are separated 
whenever possible to better isolate the frequencies of individual 
modes . 

Some definitions are helpful before proceeding with this analy- 
sis. Six subscripts which will be used are: 

S symmetric 
A asymmetric 
R right 
L left 

a “V actual or true 
m measured 

The basic symmetric and asymmetric transfer functions are: 


A YcCs) y^Cs) 


( 2 . 12 ) 


where 


yc(s) ^ i CYd C s)+y, (s)) , measurement of symmetric 
^ K L response (2.13) 


^ 7 C/j^Cs) -yj^(s)) , measurement of asymmetric 


response 


C2.14) 
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Table 2.2 

Transfer Functions of Coupled Wing Responses 
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(U|^(s)+Uj^(s) j , symmetric control input (2. IS) 

(Uj^(s) -Uj^(s) ) , asymmetric control input (2.16) 

Also, right and left control calibration errors and respec- 
tively) and the right and left measurement calibration errors (rij^ 
and rij^, respectively) are defined such that: 

actual * commanded or measured (2.17) 

^^R^measured * ^^R^ actual (2-18) 

and the left variables are related similarly. 

An example of the first case, that of exciting one wing and 
measuring one wing where no calibration errors are present, is the 
following: 



/rCs) 

= XsCs) + y^^(s) 



= 7 Tg(s) [Uj^(s)+Uj^(s)l + y T^(s) [Uj^Cs)-u^(s)] 

(2.19) 

Ur(s) 

• 1 tTs(s)*T^(s)l 

(2.20) 

/rCs) 

Hprr 

- ^ [Tj(s) - Tj^(s)l 

(2.21) 


Next, consider the last case, that of exciting both wings 
either symmetrically or asymmetrically and measuring the responses 
of both wings in the presence of both control and measurement cal- 
ibration errors. For symmetric excitation: 
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( 2 . 22 )' 


and similarly for y, and y. . From the above, it is easily shown 
that ™ ^ 


^m 


(2.23) 


and, similarly for asymmetric excitation, 


i.,. (..ip.)(..vi).,,(:!^)(«.: 


(2.24) 


The results for the remaining cases are found in a similar manner 
and presented in Table 2.2 Note that many cases follow directly 
from the above two equations when the calibration errors are set 
to zero as appropriate. For example, the case of measuring both 
wings and exciting both wings with control (but not measurement) 
calibration errors follows simply by setting rij^ = = 0 in the 

above equation. 

The following observations are made concerning these results: 


(1) Exciting one wing and measuring one wing does not allow 
the separation of symmetric and asymmetric responses. 
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(2) Exciting both wings with calibration error and measuring 
one wing means that the measured response will consist 
o£ both symmetric and asymmetric modes when only one or 
the other was wanted. A similar case exists when one wing 
is excited and both are measured in the presence o£ 
measurement error. 

(3) When both wings are excited and both measured and both 
measurement and control errors are present, then again 
the unwanted modes contribute to the response. However, 
the £ractional amount o£ this unwanted contribution is 
reduced to the order o£ the square o£ the calibration errors 

How much the presence o£ these unwanted modes a££ect the esti- 
mate o£ the damping factors of individual modes as computed from 
test data will be the discussed in Chapter III. 

2.4 MODELING OF GUST EFFECTS 

In general, there are two principal sources of random effects 
which can degrade the information content of data. These are 
process noise (e.g., gusts) and measurement noise. For a well 
instrumented aircraft, with prefiltering on data, the process noise 
effects are of most concern. Particularly for the wind tunnel tests 
at high speeds, these gust effects can obscure essential stability 
information. 

This study required a method for emulating gust effects in 
both the frequency and the time domain. Gust spectra charact is- 
tics of the 40- by 80-foot wind tunnel are largely unquantified, 
and it is not clear as to the relation between tunnel randomness 
and corresponding flight gusts. 
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The approach used is based on use of an atmospheric gust 
spectrum. The von Karman spectrum is known to be one of the most 
accurate of isotropic atmospheric gust models, although such a 
spectrum is not convenient for matching with linear filters (due 
to noninteger factors in the spectrum) . An approximating spectrum 
is the Dryden spectrum, which is "close” to the von Karman spectrum 
for low frequencies. The Dryden spectra used for longitudinal, 
lateral, and vertical gusts are: 
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where L is the correlation distance, u the wind velocity, and 

T 1 ^ 

— » T « = (w is the bandwidth of the noise). 

u^ corr '• n 

0 n 

For this work, the bandwidth of the noise was estimated at 
2 Hz (which corresponds to a correlation distance of 25.5 ft) for 
all three directions. The variance of the lateral and vertical gust 
was chosen to be a fixed fraction of the longitudinal gust variance, 
but this ratio was chosen conservatively. This allowed parameteri- 
zation of the noise-to-signal ratio (see Appendix B.2.3) as a func- 
tion of input amplitude on the longitudinal gust rms velocity. 

The resulting variance is 
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- - 2 « 

^long 


’0.485 “ 

“Lt 

- 

0.172 

-°^vert - 


_0.343 _ 


where is the total rms gust velocity. 

As discussed in Section B.2.3, these spectra and variances 
were used with the aircraft dynamic response to random gusts to 
determine the power due to random gusts in a particular measurement. 
Viewing the three gust directions as uncorrelated with each other, 
their power contributions to a particular measurement can be added 
in an rms sense. Representations of these gust effects in time 
domain simulation were approximated using white noise passed through 
a first order filter with a break frequency at 2 Hz. 

2.5 INSTRUMENTATION MODELING 

Measurements were taken from the position states except those 
measurements with which the XV-15 is now instrumented* -wing vertical 
bending acceleration and wing chordwise bending acceleration. The 
measurements were created from the appropriate linear combinations 
of the states and controls. The simulation did not include a lag 
in measurement (the strain gauge accelerometers on the pylon have 
bandwidth of 100 Hz (13.1 per rev), which is far above any modes 
of interest). Initially, random measurement noise was added to 
simulations, but in general the process noise (tunnel gust) was 
the only noise considered, due to the expected severity of such 
tunnel disturbance. 

The transfer functions all included a first order lag at 100 Hz 
(13.1 per rev), but its effect is difficult to isolate since the 
primary modes of interest lie in the regime 0.1 to 10 per rev. 
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2.6 CONTROL MODELING 


Two cases for the order of the control system actuator were 
chosen, a first order and a fourth order lag, to represent an opti- 
mistic case and worst case for the steepness with which the actuator 
gain falls off with frequency. The break frequency was 5 Ht. 

(This corresponds to a good quality control servo.) 

The frequency responses discussed in Chapters II and III all 
include a first order lag at 5 Hz. The fourth order control lag 
was achieved by applying the following third order filter to the 
transfer function data: 

ij m — with s ■ jijj. 

(s+Un) (s^+2(0.5)o)j^s+w^) 

(5 Hz = 0.654 per rev). As will be seen in Chapters III and IV, 
there is little difference in these two representations of the 
actuator below 5 Hz, for example in studying the q^^^ mode (symmetric) 
motion at 0.398 per rev). Of course, the fourth order lagged con- 
trol would make the study of high, frequency modes (greater than 
1.5 per rev) prohibitive. Cyclic inputs were studied with a first 
order control lag ac 2 Hz and 10 Hz. 

Time domain simulations of the measured responses were made 
using different inputs (multiple sinusoids, random inputs, and 
swept sine), all using a first order lag, primarily at 5 Hz band- 
width but also at 10 Hz. 

Scale errors between the commanded and actual control were 
modeled as described in Section 2.3. 




III. ANALYSIS OF ROTOR/ CANT I LEVER WING DYNAMICS 

This section presents the analytical evaluation of test con- 
siderations for rotor/cantilever wing modes. This evaluation is 
based on the frequency response characteristics of the nine degree- 
of-freedom model, discussed in Section 3.1. Section 3.2 reviews 
the wing modes which are of principal significance in defining 
stability tests at maximum tunnel speed of 190 knots. Section 3.3 
shows the input considerations which are required to isolate these 
principal modes, and Section 3.4 discusses the effect of tunnel in- 
duced random disturbances in the measured responses. A summary is 
presented in Section 3.5. 

3.1 FREQUENCY RESPONSE OF THE ROTOR/ CANT I LEVER WING MODEL 
As discussed in Chapter II, the transfer function 

^ . H(sI-F)-1g*D 

can be used to evaluate the frequency response characteristics of 
the system with dynamics matrix, F; measurement distribution matrix, 
H; control distribution matrix, G; and direct measured inputs, D. 

The stability of this transfer function is completely described by 
the roots of the characteristic equation of F. Measurement of these 
stability characteristics, however, depends not only on these roots, 
but also on the roots of the numerator of the transfer function. 
These numerator roots are governed not only by the system dynamics 
(F) , but also the mesurements and the controls which define the 
test configuration. 
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3.1.1 Transfer Functions of Rotor/Cantilever Wing Model 

The primary controls available on the XV- IS are collective 
pitch (9^) and wing flaperon (<5£), at least for the fully converted 
configuration. For reference, the frequency response to these 
inputs for all degrees of freedom of the basic rotor/cantilever 
wing, are given in Appendix C. 

3.1.2 E valuation Method of Transfer Functions 

■ I ■!■■ ■ II ■— l» ■ I ■■■■I I ■ I I II I ■ ■ I I t 

The transfer functions of Appendix C serve as a useful refer- 
ence description of the frequency response of the nine degrees -of - 
freedom to collective and flaperon inputs. In this section, we 
examine the total information content of the transfer functions, 
assuming that all states are measured and input bandwidth is above 
the highest frequency mode of the system. This "ideal" case demon- 
strates the relative ranking of modal information independent of 
test input or instrumentation limits . In Sections 3. 2-3. 4, we dis- 
cuss the effect of these limits in more detail. The purpose of 
presenting these "ideal" transfer function characteristics is to 
formulate a basis for discussing the desirability of more stringent 
requirements on instrumentation and excitation hardware. 

The basic objective of stability testing is the determination 
of frequency and damping of system responses. For a multivariable 
rotorcraft system, however, it is necessary to isolate the frequency 
and damping of the elemental modes which produce that system 
response. The problem is, however, that the system response (as 
indicated by the transfer functions of Appendix C) represent the 
sum of contributions of all these elemental modes at a particular 
frequency. In many cases, the elemental modes are essentially un- 
coupled, and a particular peak corresponds to a unique degree of 
freedom (e.g., wing vertical bending). For the rotorcraft, however, 
significant modes are highly coupled and measurement of system 
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frequency response may not provide the desired stability character- 
istics for a particular degree-of-freedom, A corollary to this is 
that determination of system damping for a particular peak, itself 
a nontrivial calculation, may not define the stability of a particu- 
lar degree of freedom. 

There are several analytical techniques which may be used to 
decompose system response into contributions from various modes 
(such decomposition is equivalent to quantifying the identif iability 
of a particular degree of freedom). These techniques include the 
following: 

(a) Mode Shape Analysis : Corresponding to a particular natur- 

al frequency of response is a mode shape associated with 
that frequency. This mode shape, or eigenvector, may be 
calculated. The eigenvector consists of a vector of 
components of the elemental degrees-of-freedom. The rela- 
tive size of these components quantifies the participation 
of each degree of freedom at a natural frequency. Unfor- 
tunately, highly coupled modes yield eigenvectors which 
show several components of nearly equal contribution, and 
it is not clear how to isolate the most significant 
degrees -of-freedom. 

(b) Residue Analysis (Appendix B) : The natural response 

of a system mode may be written as a sum of elemental 
modes. The terms of this expansion are of the form 

e ”i cos (wn^ t 

where A^^ is the combined residue of the mode eigenvalue, 

^ C^i ■ where R^^ is residue of Xj^). .Modal 

content can be estimated by ranking the residues of all 
modes at each frequency. 
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Cc) Modal Power (Appendix 3} : The modal power of a par- 

ticular transfer function peak is the power contributed 
by the complex pole pair associated with the vibration 
(in a selected frequency range about the peak) . 

(d) Modal Power Ratio (MPR) (Appendix B) : The ratio of a 

particular degrees of freedom power to total power of a 
response. The MPR can be positive or negative depending 
on whether a particular degree of freedom is adding or 
subtracting power from a resonant peak. The sum of the 
MPR* s for each mode is unity. 

The modal power and modal power ratios are used to quantify 
the following: 

(1) Which transfer function is more desirable to identify a 
particular mode? 

(2) Which mode is most significant in a resonant peak when 
more than one mode is contributing? 

The modal power ratio does not include information on the magnitude 
of the peak. (For example, transfer function B may show a resonant 
peak of much less magnitude than transfer function A, yet their 
MPR's are nearly equal.) Thus, for selecting which tran.‘:fer func- 
tion is more desirable to identify a particular mode, the modal 
power should be used. This can be done since: (1) the system ^f 
equations is normalized so the transfer function evaluations a^^e in- 
dependent of the units, and (2) the same integration interval was 
used on each transfer function. When investigating one peak of a 
particular transfer function to determine which mode is being iden- 
tified, either the MPR or modal power may be used. When comparing 
the uniqueness (meaning lack of other modes that contaminate) of 
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different peaks of one transfer function, the MPR must be used 
because the modal power may be computed at different frequency 
bandwidths for each peak. 

3.1.3 Transfer Function Comparisons 

Table 3.1 is a qualitative comparison of the transfer functions 
to collettive and flaperon. Such a table is useful for rapid refer- 
ence to determine the importance of measuring a particular degree- 
of-freedom or using a particular input. For example, the table 
indicates that the q^j^, q^ 2 » P> 2 modes can be isolated (if 
sufficient control bandwidth is available). 

Table 3.2 is the quantitative evaluation of the transfer func- 
tions, upon which Table 3.1 is based, and summarizes the modal power 
computations for the cantilever wing (symmetric motion) transfer 
functions. Transfer functions to all nine degrees -of -freedom excited 
by collective pitch and wing flaperon, as well as selected cyclic 
transfer functions, are included. In the horizontal direction of 
the table are the different transfer functions; in the vertical 
direction of the table are the different modes. As developed in 
Appendix B.2.2 and B.2.3, every mode contributes either positively 
or negatively by some great or small amount to the total power 
of each peak in a particular transfer function. Only contributions 
(in terms of MPR) greater than 51 are included in Table B.l. 

Each box associated with a particular transfer function 
contains three values. The first is the frequency of the peak 
for which the modal power is computed, and the following two values 
are the modal power and modal power ratio (MPR) . Each transfer 
function does not exhibit peaks at each of the modal frequencies. 

When a mode is not excited in a particular transfer function, the 
block for that mode is used to indicate that mode's contribution 
to the nearest adjacent peak, which is at the indicated frequency. 
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Table 3.1 

Qualitative Comparison of Transfer Functions 
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Table 3.2a 

XV- 15 Rotor/Cantilever Wing Modal Decomposition 
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Table 3.2b 

XV-15 Rotor/Cantilever Wing Modal Decomposition 
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Table 3.2c 

XV-15 Rotor/Cantilever Wing Modal Decomposition 
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Thus, the table can be used in the following ways: 

• Comparisons can be made horizontally of the modal power 
values for a desired mode to choose the transfer function 
giving the best information for that mode. This compari- 
son can be made because the equations were normalized and 
the same frequency bandwidths were used. 

• Comparisons of MPR can be made vertically to determine 
which peaks of a particular transfer function arise mostly 
from the resonance of a single mode and are least obscured 
by other modes. 

Figures 3.1a and 3.1b are illustrative of the manner in which 
modal power ratio is used to evaluate system response. The wing 
chordwise acceleration to wing flaperon transfer function 

is shown in Figure C.2 (Appendix C) . Four distinct peaks are evi- 
dent. Figure 3.1a shows the modal power ratios for each of t.iese 
peaks for contributors above 5% of the total power in the response. 
(Note that other degrees of freedom not shown are below 5% of the 
total power and sura with the plotted MPR’s to unity.) The following 
Figure 3.1b shows the effect of additional control system lag beyond 
5 Hz and the alterations in MPR associated with this lag. 

The modal power ratio has even more utility when considering 
the coupled wing responses. This will be discussed in Chapter IV, 


3.2 DISCUSSION OF TRANSFER FUNCTION CHARACTERISTICS 


The primary modes of interest for XV-15 structural stability 
testing are those associated with the wing. For the purposes of 
this study, these modes include the wing vertical bending (q^j^) , 
wing chordwise bending wing torsion (p) modes 


These 
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Figure 3.1 Modal Power Ratio; Wing Flaperon to Wing Chordwise Acceleration 
Cantilever Wing Transfer Function (Symmetric Motion) 
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modes are lightly damped and, therefore, of significant interest 
for flight test prediction. 

In addition to the wing modes, there is another lightly 
damped mode whose response is important. This is the upper inplane 
mode, although this mode is at a high frequency. These four 

modes are summarized with respect to frequency and damping in 
Table 3.3. The data in Table 3.3 is the "ideal" frequency and 
damping, based on an analysis of the characteristic equation. Sta- 
bility parameters determined from various transfer functions will, 
in general, be different, as discussed below. 

Of these four modes, further analysis was focused primarily 
on the q^j^ and q ^2 niodes. The p mode has the greatest damping of the 
four at 0.058, and the analysis performed on the q^^^ and q ^2 niodes 


Table 3.3 

Principal Symmetric Modes 


MODE 




2.43 Per Rev 
(18.6 Hz) 

0.02990 

p 

1.34 Per Rev 
(10.3 Hz) 

0.05783 

%2 

0,666 Per Rev 
(5.09 Hz) 

0.04258 

%1 

0.398 Per Rev 
(3.04 Hz) 

0.04798 
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can be easily extended to the higher two modes when consideration of 
control excitation at their frequencies is desired (p at 10.3 Hz 
and 18.6 Hz) . 

Observations about the transfer functions are made here pri- 
marily to point out modes not available in the wing bending transfer 
functions from instrumentation presently planned for the XV- 15 
(q^l/6f and presently available). For example, it is 

difficult to obtain information on the upper inplane mode. Direct 
measurements of the inplane motion show a resonant peak in the fre- 
quency vicinity of the upper inplane mode (^^^) . However, the rotor 
coning mode CB) is very close to this frequency. Excitation of 
the inplane degrees of freedom by collective gives a response with 
magnitude lO’^, but the collective also excites the 3 mode consid- 
erably (see Table 3.2). Excitation of the inplane degrees of freedom 
by the flaperon does not excite the 6 mode significantly; however, 
the resonant peak has a magnitude of only 10 . Cyclic inputs would 

excite the mode sufficiently without significant participation 
of the 8 mode. Resonant peaks with good information on the 
mode occur in the transfer functions of wing torsion (p) and rotor 
inplane motion and B^g^) to cyclic inputs and ‘ 

Thus, by measuring wing torsion, it would not be necessary to add 
strain gauges on the rotor blades to measure the mode. 

Collective pitch to the rotor coning degree of freedom does 
not show sharply distinguished resonances for modal decomposition. 
Excitation of this degree of freedom with wing flaperon does show 
resonance; however, the mode corrupts the considerably (see 
Table 3.2). Also, this peak is rather small in magnitude (10 ^) , 
again indicating that a rotor blade measurement would not add to 
the information available from present instrumentation. 
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Gimbal measurements excited by flaperon (6£) give no additional 
modes. Collective pitch (9^) does excite the mode: however, it 
is corrupted somewhat by the B mode (see Table 3.2). The rotor 
azimuth perturbation transfer functions give no additional modes. 

Of the wing bending measurements presently instrumented on the 
XV- 15, investigation of the cross transfer functions shows that 
q^z/'Sf quite attenuated, with peaks at a magnitude of 10*^ or 
below. However, actually more desirable than 

for identifying the q^j^ mode, showing a greater magnitude, greater 
modal power and also a greater modal power ratio (Table 3.2). Both 
measurements, q^j^ and q^2’ recorded when either 9^ or 

is excited. Obtaining the cross transfer functions (q^j^ /e or 
q^2/'5r) is merely a matter of additional data processing. For iden- 
tifying the damping of the q^^ mode, the estimates from the 
and transfer functions can be combined to give an improved 

estimate. The transfer function can be used to identify 

the p mode. (There are other transfer functions where the p mode is 
excited more such as the gimbal angle, inplane bending, and torsion 
transfer functions^ but these require additional instrumentation.) 
The function can be used to identify the 6 mode 

(given that the control system could be excited at this high fre- 
quency: 19.5 Hz). As with the p mode, the B mode could be identi- 
fied from other trar 'fer functions fay adding measurements. The 
wing bending measurements excited with cyclic yield no additional 
modes + and 8 mode obscure each other). 

Torsion measurements excited by either collective pitch or wing 
flaperon give no additional modes. The significant observation 
about a torsion measurement, however, is that the upper inplane 
mode can be identified with excitation by cyclic pitch (if the 
control system could be excited at this frequency: 18.6 Hz). .A 
torsion measurement could be obtained from linear combinations of 
vertically aligned ar.celerometers on the fore and aft ends of the 
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pylon. This instrumentation would be much easier than adding a 
strain gauge to the rotor blades and passing the signal through 
a commutator at the rotor hub. 

3.3 INPUT FREQUENCY REQUIREMENT 

The overall system transfer function analysis of Sections 3.1 
and 3.2 is now specialized to analysis of the wing transfer func- 
tions. From the results of Tables 3.1 and 3.2, wing response trans 
fer functions are isolated which would give the best information 
about particular degrees of freedom. 

The damping of modes from particular transfer functions is now 
presented from the frequency response data plotted in Figures C.l 
to C.12. The method used is that suggested by Johnson in Ref. 6, 
Appendi.x E. A detailed derivation and explanation of the method is 
included in Appendix B of this report. 

Table 3.4 shows the damping factors calculated for each 
mode using this method. The transfer function was selected using 
the modal power as a criterion. The table includes the effects of 
two different control lags- -a first order and a fourth order lag-- 
both with a bandwidth of 5 Hz. It is seen that, with existing 
wing input channels, wide band input frequencies can be used to 
isolate coning (8^), upper blade out-of-plane and inplane 

modes, and wing torsion. Only wing vertical bending and 
chordwise bending can be excited with a one/rev limit on 9^ and 5^. 
Note that the mode could be excited by a cyclic input (although 
at a high input frequency). 

The calculation of the dampings of Table 3.4 was achieved by 
integrating over the frequency band 0.9 to 1.1 where is 

the peak frequency. Damping obtained by this method was found very 
sensitive to this frequency band. A wider band, 0.8 • 1.2 

was found too wide, introducing other modes and degrading the damp- 
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ing to S0*<! of the theoretica'' <'lue for the and modes. 
Reducing the frequency bandwi ' ‘'.J' to the half -power value of 0.9S 
1.05 further showed a degr'vjd accuracy. 

The response conclusions based on the transfer functions apply 
principally to input and measurement frequency requirements. Ampli- 
tude specifications for inputs and instrumentation are based on 
actuator and sensor limits not fully available at this time, but 
easily determined from the transfer functions of Appendi.x C. An 
important additional source of input requirement arises from the 
consideration of tunnel induced turbulence, discussed in the follow- 
ing Section 3.4. 

3.4 THE EFFECT OF NOISE IN SELECTING INPUT AMPLITUDES 

Having selected the transfer functions to identify the primary 
modes on a deterministic basis, we now address the effect of random 
disturbances in the wind tunnel in further specifying input ampli- 
tudes. There are two questions of interest here: 

(1) Ifhat is the responsiveness (or susceptibility to dis- 
turbance) of a particular measurement to random gusts? 

(2) Given a particular measurement e.xcited by a particular 
input, what input magnitude is required to achieve a 
desired signal-to-noise ratio (SNR), or, alternately, 
minimize a noise-to-signal ratio (NSR)? 

These questions can be answered by considering the model of 
the XV-IS dynamics 
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X • Fx + Gu + Tv 
^ ■ Hx + Du 

We have been considering the steady state response to control exci- 
tation 


y(ju)) « [HCju)I-F)'^G + D] u(jo)) 

The resulting Bode plot can be used to predict the magnitude and 
phase of the response to sinusoidal excitation (a discrete fre- 
quency) or band limited random excitation in a frequency range of 
interest. For example, the entire Bode plot would be the response 
to excitation of all frequencies in the range of the Bode plot or 
white noise limited to this range. Thus, we can consider random 
excitation in any frequency band; and, specifically, we can consider 
random excitation in the vicinity of a peak of interest. 


It may be shown (Appendix B.2.3), that the noise-to-signal 
ratio in a frequency band ^2 - can be written 


N 

5 


r(u)2-wi) 

uTT ^ 


T 


rms s rms s 


where 


(3.1) 


is the noise power in the frequency band 

Xg is the equivalent power from a sinusoidal signal in 
the same frequency band 

q is the gust power spectral density 

r is the measurement noise power spectral density 

^rms root -mean- square value of the input 
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The second terra of this equation is negligible corapared to 
high speed wind tunnel turbulence, and may be neglected. (Note 
that the effect of measurement noise can be simply added.) 

Plots for the selected transfer functions of NSR, as a func- 
tion of the rms input amplitude, are plotted as families of curves 
parameterized on tue rms velocity of the gust in the wind tunnel 
and on the control actuator bandwidth (Figures 3.2 through 3.7). 

The results are presented for two cases. 

(1) First order lagged control with the Dryden wind gust 
model (see Section 2.4). 

(2) Fourth order lagged control with the Dryden wind gust 
model . 

Two general observations are evident in the NSR plots (Figures 
3. 2- 3. 7). First, the 2 Hz control bandwidth case shows a greater 
deviation from the base case of 5 Hz than the 10 Hz control band- 
width case. This is to be expected since the £ Hz break frequency 
is close to the upper limit of the wing modes under consideration, 
while the 2 Hz break frequency is completely below them. Second, 
the fourth order lagged control system makes the break frequency 
more significant, and indicates a significant increase in rms 
amplitude to achieve a desired NSR for a mode above the break fre- 
quency. 

An example use of these charts is as follows. To achieve a 
NSR between 0.1 and 0.2 for the q^2^®o ^^®^sfer function when the 
expected total rms wind tunnel gust is 10 ft/sec, the following 
collective pitch input would be required: 
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Tilt Rotor at 190 knots 
First Order Control Lag at 5 Hz 
First Order Measurement Lag at 100 Hz 
Oryden Gust Model 
Correlation Distance 25.5 ft 


RMS Gust Velocity In 


CONTROL BANDWIDTH: 

2 Hertz 

5 Hertz 

10 Hertz 


Figure 3.2 Noise-to-Signal Ratio- -Collective Pitch to Wing Chord 
wise Acceleration Cantilever Wing Transfer 

Function (For Frequency Range 4.20 to 6.35 Hz) 
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Tilt Rotor at 190 knots 
Fourth Order Control Lag at 5 Hz 
First Order Measurement Lag at 100 Hz 
Dryden Gust Model 
Correlation Distance 25.5 ft 



Figure 3.3 Noise- to-Signal Ratio- -Collective Pitch to Wing Chord- 
wise Acceleration Cantilever Wing Transfer 

Function (For Frequency Range 4.20 to 6.35 Ht) 
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RHS Gust Velocity In 


Tilt Rotor at 190 Kts 
Fourth Order Control Lag at 5 Hz 
First Order Measurement Lag at 100 Hz 
Oryden Gust Model 
Correlation Distance 2S.S ft 


CONTROL BANDWIDTH: 

2 Hertz 

5 Hertz — 

10 Hertz 


if, DEGREES 


Figure 3.4 Noise-to-Signal Ratio- -Wing Flaperon to Wing Vertical 
Acceleration Cantilever Wing Transfer Func- 

tion (For Frequency Range 3.04 to 4.60 Hz) 
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Y 10 \l0— RMS Gust Velocity in 
\\ \ ft/sec 


Tilt Rotor at 190 knots 
First Order Control Lag at 5 Hz 
First Order Measurement Lag at IOC Hz 
Oryden Gust Model 
Correlation Distance 25.5 ft 


CONTROL BANDWIDTH: 

2 Hertz 

5 Hertz 

10 Hertz 


if, DEGREES 


Figure 3.5 Noise- to-Signal Ratio--Wing Flaperon to Wing Vertical 
Acceleration Cantilever Wing Transfer Func- 

tion (For Frequency Range 3.04 to 4.60 Hz) 




Tilt Rotor at 190 knots 
First Order Control Lag at 5 Hz 
First Order Measurement Lag 100 Hz 
Oryden Oust Model 
Correlation Distance 2S.S ft 


RMS Gust Velocity In 
/sec 


CONTROL BANDWIDTH: 

2 Hertz 

5 Hertz 

10 Hertz 


DEGREES 


Figure 3.6 Noise-to-Signal Ratio- -Collective Pitch to Wing Verti 
cal Acceleration Transfer Function for Asym 

metric Motion (For Frequency Range 3.33 to 5.29 Hz) 
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CONTROL BANDWIDTH 
(HERTZ) 

INPUT AMPLITUDE 
IN DEGREES - RMS 
(FIRST ORDER UG IN CONTROL) 

2 

5.80 - 8.10 

5 

2.97 - 4.20 

10 

2.31 - 3.30 


(From Figure 3.2) 


The NSR plots can be used in the following manner (neglecting 
measurement noise) : 

(1) A range of acceptable SNR’s and, hence, NSR's, and a 
range of expected rms longitudinal gust should be 
selected. 

(2) The corresponding limiting upper and lower rms input 
amplitudes from the NSR plots give the required ampli- 
tude to achieve the desired SNR. 

(3) These values should be determined for the first order 
lagged control model and the fourth order lagged model 
(at the desired break frequency) to give a best and worst 
case (as discussed in Section 2.6). 

3 . 5 SUMMARY 

This section has presented the results of a detailed study of 
a model of XV-IS rotor/cantilever wing dynamics at 190 knots. 

The objectives of this study have been the following: 

(1) Calculation of transfer functions from the two principal 
controls (collective and flaperon) to the nine degre-^s- 
of-freedom of the model. 
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(2) Ranking of the transfer functions with respect to their 
utility in isolating principal modes, and input frequency 
bandwidth requirements. 

(3) Quantification of the effects of tunnel noise on required 
input amplitudes. 

It is concluded that the principal limitation to mode isolation, 
for the nine degree-of-freedom model is input bandwidth. Instru- 
mentation of the wing is sufficient to identify all but the 
and lower (-1) rotor modes if collective and flaperon input channels 
are used. To identify the mode, measurements of wing torsion 
excited by longitudinal cyclic would be effective. In all cases, 
however, input frequencies required to identify the upper (+1) 
rotor modes and wing torsion mode are above 1 per rev. The lower 
(-1) rotor modes are not identifiable with wing instrumentation 
because they are heavily damped and obscured by other modes. 

Evaluations of the quality of transfer functions from this 
nine degree-of-freedom model are based on a control system band- 
width of 5 Hz. A first order and fourth order lag were used with 
this bandwidth. The principal effect of increased control system 
rolloff beyond 5 Hz is increase in input a.i'.plitude required to 
achieve an adequate noise- to-signal ratio of output data. Input 
amplitudes required to achieve a desired noise- to-signal ratio 
for different levels of wind gust were determined for a control 
system bandwidth of 2, 5, and 10 Hz. 
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IV. COUPLED WING DYNAMICS 


This section presents the analytical evaluation of the coupled 
wing dynamics. Section 4.1 discusses the frequency response for 
the measurements presently available (both wing bending accelera- 
tions) excited by collective pitch and wing flaperon. Section 4.2 
summarizes the modeling of calibration errors in measurement and 
excitation and discusses the results. Section 4.3 presents the 
effect of tunnel -induced random disturbances and the input ampli- 
tudes and frequency ranges whicn are required to identify principal 
modes of interest with particular noise-to-signal ratios. 

4.1 FREQUENCY RESPONSE OF THE COUPLED WING MODEL 


4.1.1 Computation of Coupled Wing Responses 


The discussion of the coupled wing model will be limited to 
the and Q;^2^®o ^^^^sfer functions for identification of the 

following modes: 


q^^g! symmetric vertical wing bending mode 
q,^: asymmetric vertical wing bending mode 

q 2 g: symmetric chordwise bending mode 

q 2 y^: asymmetric chordwise bending mode 


As discussed in Section 2.3, the response of the right wing by 
excitation of the right wing has a symmetric and an asymmetric por- 
tion. Excitation of a wing, and measurement of that wing, is not 
sufficient to separate the symmetric and asymmetric coupled wing 
modes if they are "close in fre.; . ." As discussed in Section 

2.3 and again in Section 4.2, additional instrumentation makes it 
possible to separate the symmetric and asymmetric frequency 
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responses, except for some fraction of the opposite motion (asym- 
metric or symmetric) due to calibration errors. 

The nine degree-of -freedom rotor/cantilever wing model can 
be used to analyze the coupled wing (e.g., right and left wing 
coupled modes) in the following manner: 

(1) compute the symmetric transfer function, Tg(ju), from 
the linear equations representing the symmetric motion; 

(2) compute the asymmetric transfer function, T^Cju), simi- 
larly, and 

(3) combine these transfer functions at discrete values of 
the frequency, cu, by the formula 

T(jw) - ^ (TgCjo)) ♦-T^(jw)] . . 

Because this new transfer function T now has twice as many 
modes in the same frequency range, more pronounced modal coupling 
occurs in a small frequency range and the phase may appear to be 
discontinuous (Appendix B.2.1). However, if a sufficiently small 
frequency increment were used, it would show the phase to be con- 
tinuous, since both transfer functions are analytic. 

Clearly, the coupled wing motion will exhibit characteristics 
of both the symmetric and asymmetric motions. The characteristics 
of the symmetric transfer functions were discussed in some detail 
in the previous chapter. Before continuing the discussion of the 
coupled wing motion, the characteristics of asymmetric motion will 
be presented briefly. 
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4.1.2 Asymmetric Mode Characteristics 

The principal asymmetric modes are listed in Table 4.1. The 
asymmetric frequency responses for the wing bending accelerations 
are shown in Appendix C (Figures C.2S to C.28). 


The linear model for the asymmetric motion differs from the 
symmetric motions in that the coupling due to the drive shaft 
between the two rotors introduces an oscillatory mode associated 
with the rotor speed perturbation. The wing .stiffnesses, associated 
with the two asymmetric wing bending degrees of freedom, were in- 
creased to shift the and q 2 ^ modes up in frequency from the 
and q 2 g modes in order to distinguish them. The q 2 ^ mode was 
shifted up such that it was almost coincident with the mode 
(the rotor speed perturbation mode). This was considered one type 
of worst case for analyzing the effect of the asymmetric mode on 
measured frequency and damping. 


Table 4.1 

Principal Asymmetric Modes 


MODE 

DAMPED 

FREQUENCY 

DAMPING 

FACTOR 

C 

^+1 

2.44 Per Rev 
18.6 Hz 

0.02978 

P 

1.42 Per Rev 
10.8 Hz 

0.05382 

*^w2 

0.727 Per Rev 
5.55 Hz 

0.02124 


0.562 Per Rev 
4.29 Hr 

! 

0.05543 
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4.1,3 Coupled Wing Mode Characteristics 

The transfer functions for the coupled wing responses are calcu- 
lated from the plots presented in Appendix C. For the two principal 
modes of most interest, and q^ 2 * *^1A prac- 

tically unobscured and the q 2 g and q 2 ^ modes are obscured to a noder- 
ate and greater extent, respectively. This illustrates both ex- 
tremes which can occur in the actual dynamics of the XV-IS. 

The transfer function shows a prominent resonant peak 

for q 22 (at 0.66 per rev), but the q 2 ^ mode (at 0.72 per rev) 
shows no resonant peak at all. It is obscured by the q 2 g and 
modes. The q 2 g peak itself is obscured somewhat by the mode 
(at 0.73 per rev), as is shown by the modal power ratio (MPR) in 
Figure 4.1. The symmetric and asymmetric flapping modes Bg and B^ 
are very close in frequency (Bg ■ 2.55 and B^ ■ 2.61 per rev), 
and although they are not principal modes, it is interesting to 
note their equal contributions to the second resonant peak on tl 

^^**’'sfer function from the MPR shown in Figure 4.1. Because 
the symmetric and asymmetric B modes result in a single resonant 
peak, these modes cannot be distinguished in this transfer function, 
which illustrates the value of being able to separate the symmetric 
and asymmetric motions. 

The transfer function shows distinct resonant peaks 

for the qj^g and q^^^ modes and one resonant peak where the two 
torsion modes occur (pg at 1.34 and at 1.417 per rev). The 
MPR's for this transfer function, which are shown in Figure 4.2, 
indicate that q^^g is unobscured and that q^^^ is slightly obscured 
by qj^g. The torsion modes obscure each other in the third peak 
as expected. The modal power ratio for collective pitch to wing 
vertical bending acceleration (q^^/Qg) with coupled wing motion is 
shown in Figure 4.3. 

At this point, it must be concluded that the difficulties 
with the coupled wing model are as follows: 
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Figure 4.2 V.odal Power Ratio (MPR)--Wir.g Flaperon to Wing Vertical 
Bending .4cceleration Coupled Wing Motion 










gure 4.3 Modal Power Ratio (MPR)--Wing Vertical Bending Acceler- 
ation to Collective Pitch Coupled Wing Motion 
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(1) There are so many modes in the same frequency region 
that it is difficult to determine whether a particular 
resonant peak is due to symmetric or asymmetric motion 
(i.e., due to the excitation of the wing which is being 
measured, or the reaction force of the rest of the air- 
craft on the wing being measured). 

(2) Even more importantly, given a particular resonant peak 
on the frequency response, one does not know how many, 
which one, or how much other modes are contributing to 
that peak. The modeled dynamics used in this analysis 
show examples where some principal modes are obscured 
and others are not. The measured dynamics of the XV- 15 
could be better or worse, but the significant point is 
that one will not be able to determine this unless the 
symmetric and asymmetric motions are separated. It is 
important to know whether a principal mode is obscured 
or not and to know how much. The damping calculated 
from a resonant peak could show the XV-15 to be quite 
stable, when in fact there are two or more modes contrib- 
uting to this peak, one of which is safely stable, but 
another of which has very light damping that could become 
unstable in some flight conditions. 

4.2 EVALUATION OF CALIBRATION ERRORS 

The approach for modeling calibration errors for various wii. i 
tunnel test configurations C©‘g-» excite one wing and measure one 
wing, or excite one wing and measure both wings, etc.) is developed 
in Section 2.3. In this section, throe test configurations are 
considered in detail to assess the effect of unwanted modes on 
the damping estimates' of desired modes. The three cases are: 


(1) the present test configuration (excite one wing and 
measure that wing) ; 
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(2) An improved test configuration which will separate sym- 
metric and asymmetric motion (excite one wing and 
measure both wings) ; and 

(3) the most accurate test configuration for separating sym- 
metric and asummetric motion (excite both wings and 
measure both wings). 

For each configuration» the actual environment is considered; that 
is, calibration errors in both measurements and controls are 
modeled. 


First consider the effect of calibration errors in both 
measurement and control of the present test configuration. The 
notation of Section 2.3 is used here except that functional 
dependence on the Laplace complex variable s is not shown 
explicitly. The measured response of the right wing is 


yR ■ * ’'i) 

m a 

where is the measurement calibration error and the actual 
response of the right wing is 


s “s 


Va. 


Thus , 


^R„ ■ * Va > 

m a a 

* ^A |t ■ X'^ * "l’] 
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where and are the control calibration errors, 
y R 

^ - ^ (1 > nR)Cl * eR)[Ts * T^l 

» [1 + £j^ + Hr + ej^l I [Tg + T^] 

Thus, calibration errors in both excitation and measurement for 
the present configuration have the effect of scaling the measured 
transfer function by a gain factor, which does not affect the 
determination of the resonant peak frequencies. "This gain factor 
also does not affect the damping calculation because the square 
of the gain enters in both the numerator and denominator and, there- 
fore, cancels out. 

Thus, for the present wind tunnel test configuration, calibra- 
tion errors have a negligible effect in the calculation of the 
frequency and damping associated with a particular resonant peak 
on the frequency response. However, the difficulty with the present 
test configuration, as discussed in the previous section, is that 
the symmetric and asymmetric modes obscure each other, making the 
determination of the damping of principal modes for assessing the 
dynamic stability of the XV-15 ineffective. 

Therefore, it is very desirable to separate the symmetric and 
asymmetric motion to determine with .confidence which modes are 
present in the resonant peaks of the frequency response. This can 
be accomplished in three ways, as discussed in Section 2.3: 

(1) Measure both wings exciting one wing. 

(2) Measure one wing exciting both wings. 

(3) Measure both wings excising both wings. 
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It will be shown later that the third alternative, although it in- 
volves more test equipment and, hence, more cost, actually reduces 
the effect of calibration errors to a negligible level. The first 
alternative will be considered here since it is the easiest to im- 
plement. 


Again, using the notation of Section 2.3, 


m m 


Considering the case where the control input is to the right wing 
only (i.e. , Uj^ ■ 0) 

Under the assumption that products of errors may be neglected. 


m ' 


(4.1) 


Likewise, 






z 





(4.2) 


Since a gain constant times the measured transfer function does not 
affect the ability to identify the frequency and damping of prin- 
cipal modes, Eqs. (4.1) and (4.2) can be written as: 
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(4.3) 


(4.4) 


Let 


A ^ ^R'^l 


(4.S) 


be called the "distortion factor" because it is this factor times 
the asymmetric transfer function in Eq. (4.3) which "distorts" the 
desired symmetric transfer function. Similarly, this same factor 
times Tg distorts in Eq. (4.4). That is. 



oc 


Ts * 


C4.«) 




OC 




(4.7) 


For values of the measurement and control scaling factors 
between -0.20 and +0.20, the distortion factor varies in the range 


-0.25 < ^ ± +0.25 

Clearly, the value of the distortion factor is most strongly 
dependent upon the difference of the calibration errors in the 
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right and left wing response measurements, less upon the 

right control calibration error, £j^. Figure 4.4 shows how A varies 
as a function of the right and left measurement calibration errors 
for discrete values of the control calibration error. 

To study the effect of calibration errors on the damping esti- 
mates, the mode was chosen since both the symmetric and 
asymmetric modes are not obscured significantly by other modes. 

The symmetric or asymmetric transfer function was degraded by adding 
or subtracting the other transfer function with varying levels of 
distortion factor. The dampings calculated for these various levels 
of distortion factor are plotted in Figure 4.5. As can be seen, 
calibration errors do not have a significant effect on the calcu- 
lated values of the q^^g and q^^^ modal damping. The error in the 
calculated value is primarily due to the influence of other modes 
which are present even when the calibration errors are zero. 

There are two important considerations to be kept in mind. con- 
cerning the effect of calibration error on the calculated damping 
values. First, the calculated damping values presented here are 
empirical in the sense that the complete algorithm for accurately 
determining c is a subject for further study. Naturally, there are 
several data processing considerations affecting the accuracy at c 
(see Appendix B for a discussion of the method used to calculate 
damping). These considerations include the distance between the 
frequency response data points (a function of the data record 
length), the size of the bandwidth about the peak used in the cal- 
culations, and with actual data, the signal-to-noise ratio of each 
data record and the number of data records averaged. It was found 
that the damping value (;) was somewhat sensitive to the size of 
the bandwidth chosen about the peak. Too wide a bandwidth included 
too much influence from other modes and tended to underestimate ;; 
too narrow a bandwith did not contain sufficient information and 
also underestimated Better damping values could be produced by 
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Figure 4.S The Influence of Calibration Errors on the Calculated 
Values of Damping for the q - Mode from the cj ./9 
Transfer Function ^ ^ 
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«. calculating the transfer function with smaller frequency increments, 

and varying the bandwidth used for each peak (within what is known 
to be a reasonable range from the plot of the frequency response) 
to achieve an accurate value for c. The purpose of results pre- 
sented here is to provide sufficient information to evaluate vari- 
ous test instrumentation configurations. 


A second problem in analyzing the effect of calibration errors 
on calculated damping values is knowing what mode is actually being 
measured from a particular resonant peak. The q ^2 asymmetric (Q2 a^ 
peak gives an excellent example of what happens to the damping 
calculation when a mode is being obscured by other modes. The q 2 ^ 
mode (at 0.72 per rev) is obscured somewhat by the q 2 g mode (at 
0.66 per rev) and significantly by the rotor speed perturbation, 
which is an asymmetric mode (at 0.73 per rev). It is not valuable 


to plot the damping of the q 2 ^ mode as a function of 
errors since it is obscured by other modes. However 


calibration 
the calculated 


values for the damping for different levels of distortion are given 


in Table 4.2 for reference. 


Table 4.2 

Calculated Values of Damping for q 2 ^ Mode 
from q Transfer Function 

W L 0 


DISTORTION 

CALCULATED 

FACTOR 


0.0 

0.0429 

0.25 

0.0439 

1.0 

0.0426 

Actual Value 

of ; » 0.0212 
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The effect of calibration errors on the most complete test in- 
strumentation configuration will now be discussed. Measuring both 
wings while exciting both wings with calibration errors on both 
measurement and control will be called the full test instrumentation 
configuration. 

It follows from the expression for the theoretical transfer 
function in the above case in Table 2.3.1 that 



k 


Ur-Sl) Cir -n,^) \ 

l2*nR+n[_J Aj 



k 


(sR-El_)(nR-nL) 


+ 


I 


where 

k - 1 * 7 t^R+nj^tZER) 


The distortion factor here becomes 




Considering that each calibration is 20“ and that they are combined 
in the worst case, the magnitude of A is bounded at 4%. For the 
previous test configuration (measure both wings, excite one wing) 
the magnitude of A was bounded a*. 25%. Thus, by exciting both 
wings, even though an additional source of calibration error is 
introduced, the maximum distortion factor is reduced by a factor 
of six. Based on the variation observed in the damping as a func- 
tion of the distortion factor in the previous discussion, it is 
concluded that the effect of calibration error on damping determin- 


es 


! 


ations is negligible using this full test instrumentation config- 
uration. 

4.3 REQUIRED INPUT AMPLITUDES TO ACHIEVE A DESIRED SIGNAL-TO-NOISE 

RATIO FOR THE COUPLED WING 

The input amplitude necessary to achieve a desired noise-to- 
signal ratio will be discussed using the vertical wing bending sym- 
metric and asymmetric modes as examples. It is assumed 

that the symmetric and asymmetric motions are not decoupled; that 
is, the transfer functions from which these results were obtained 
are the sum of the symmetric and asymmetric transfer functions. 

This simulates the case when only one wing is excited and only one 
wing measured. 

The noise- to-signal (NSR) ratio calculation for a particular 
mode is based upon the frequency response to the desired control 
and the frequency response to gusts in a frequency band about the 
resonant peak for that mode. This calculation is described in 
Appendix B. For example, assume information on the q^^g mode is 
sought from the frequency response of the vertical bending accelera- 
tion to flaperon input (i.e., the ^’^^’^sfer function). The 

bandwidth about the q^^g resonant peak for which the NSR calcula- 
tions apply is 2.53 to 3.83 Hz (0.335 to 0.501 per rev). 

The result is shown in Figure 4.6. From this figure, it can 
be seen, for example, that, if the root-mean- square gust velocity 
is 1 ft/sec, a noise- to- signal ratio of 0.1 would require 3.6 de- 
grees rms flaperon input. As another example, for a gust velocity 
of 10 ft/sec rms, the same flaperon input results in a noise-to- 
signal ratio of 10. (These examples are for a first order control 
lag with a break frequency at 5 Hz). 
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Noise-to-Signal Ratio for Symmetric Verticil Wing Bend 
ing (q^g) Mode from Transfer Function (For Fre- 
quency Range 2.53 to 3.83 Hi: About the q^g Peak) 







Next, assume information on the mode is sought from the 
same transfer function. Therefore, the calculations are 

based on a bandwidth about the resonant peak (3.83 to 4.82 Hz) 
The results are shown in Figure 4.7. 


A comparison of Figures 4.6 and 4.7 shows that the same rms 
flapercn input in the presence of the same rms gust velocity would 
mean a higher (i e. , worse) noise-to-signal ratio (.NSR) for the q^^ 
mode than for the q,^ mode. This was to be expected because the 
q^^^ resonant peak has a lower amplitude than the qj^g resonant peak. 
Large inputs are required to lift the signal above the noise. 

A more favorable situation e.'r'.sts when information on the q^^ 
mode :s sough; from the response of q.., to collective pitch, inputs 
because the resonant peak of the mode has a higher amplitude 
in the transfer function. The NSR results for this case, 

shown in Figure 4.8, are considerably improved over the results in 
Figure 4.7, which were based on the transfer function. 


4 . 4 SUMMARY 

The test configuration of exciting one wing and measurrng 
that wing does not allow separation of the symmetric and asymmetric 
modal responses. It follows t’ t determination of frequer<'v and 
damping of particular modes is difficult, requiring very special- 
ized input designs or more complicated instrumentation. 

If both wings are exci-ced and/or the res> onses of both wings 
are measured, symmetric and asymmetric motions can be separated. 

In this configuration, control and measurement calibration errors 
will result in some fract.onal part of unwa.^ted modes tc be present 
with desired modes. 
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Figure 4.7 Noise-to-J: ignal Ratio lor Asymmetric Vertical Wing Bend 
ing Mode from Transfer Function (For Fre- 

quency Range 3.83 to 4.82 Hz About the Peak) 
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Figure 4.8 Noise- to-Signal Ratio for Asymmetric Vertical Wing Bend 
ing Mode from Transfer Function (For Fre- 

quency Range 3.49 to S.29 Hz About the q^^ Psak) 
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Measurement and control calibration errors less than 20% do not 
produce significant error in calculating damping values compared 
to the attenuation effects of wing cross-coupling. 
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V. EXAMPLES OF EV/LUATION OF METHODS FOR 
MODAL IDENTIFICATION FROM DATA 


The preceding two sections have summarised analytical results 
on input and measurement system requirements for XV- 15 wind tunnel 
tests. For any specified test configu:c?tion , however, the final 
accuracy of the modal frequency and damping from test data depends , 
on the processing methods used. This section discusses a prelim- 
inary evaluation of two such methods - -spectral analysis and maximum 
likelihood parameter identification. 

The objective of this phase of the study was twofold. First, 
it was desired to use the digital simulation as a data generator 
for evaluating algorithms to obtain more accurate estimates of 
frequency and damping. Second, use of an advanced parameter iden- 
tification algorithm was to be performed on such data, and compared 
with the more conventional spectral analysis approaches. 

5.1 SUMMARY OF ALGORITHMS 


The spectral analysis of this data was performed by a Time 
Series Analysis Program, which computes an estimate of the transfer 
function, H(f) , by dividing the estimated cross-spectrum of the 
input and output channels by the estimated auto-spectrum of the 
input channel. That is, 


H(f) 




The spectra were computed by the Fast Fourier Transform algorithm. 
Details of this computer program are found in Reference 7 . 


1 
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The maximum likelihood parameter identificat • on was performed 
with a modification of an existing program, SC^DNI [8]. The modi 
fication was to use a simplified model of a second order frequency 
damping model with measurements only of acceleration. Note that, 
in general, it is necessary to perform a determination of transfer 
function order prior to actual identification. Development and 
programming of an analysis to do this for the current problem is 
beyond the scope of this program, and so the assumption of the 
second order model was used. For the purposes of this study, this 
model was considered satisfactory. 

The model was as follows : 


^l' 


■ 0 

0) 




"0* 


’Xj^CO)' 

m 

S5 





+ 


u, x(0) = 


— 1 

•X 

1 


-w 

-2?oj 


.^ 2 . 


_g_ 


_X2(0)_ 



where the parameters to be identified are w, g, and possibly 
the initial conditions, Xj^CO) and X 2 ( 0 )* 

Note that this model differs from the standard state variable 

2 

model for a second order system (where F ^2 * ^ ^21 * ”‘** ) • 

form is more amenable to identification than the standard canonical 

form. 

5.2 SIMULATED TEST DATA 

The data were generated by a digital computer simulation based 
on the XV-15 mathematical modex discussed in Chapter II. It was 
assumed that a means existed which allowed separation of the sym- 
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metric and asymmetric responses of the coupled wing motion. That 
is, either both wings were excited or the responses of both wings 
were measured. Hence, the simulation used the equations of motion 
for the nine degree-of-freedom symmetric model. Two test cases 
were generated. In the first simulation case, a Gaussian random 
input in collective pitch was passed through a first order lag with 
a break frequency of 5 Hz (0.654 per rev). This lag was introduced 
to approximate the dynamics of the control system. The data length, 
was 20 revolutions of the rotor, and data samples were taken every 
0.04 of a revolution (0.25 rad) r 

In the second case, the collective pitch input was the sum 
of five sine waves whose frequencies were in the neighborhood of 
the resonant peak (specifically, they were 0.615, 0.628, 0.698, 
0.739, and 0.785 per rev). The steady state response to this 
input was simulated for 12 revolutions of the rotor and the data 
sampled every 0.0239 of a revolution (0.15 rad). 

5.3 RESULTS 

The spectral analysis method was only used on the first set 
of test data (e.g., random input). When applied to the measurement 
of wing chordwise acceleration, an estimate of the ^^^i^sfer 

function was obtained. It displayed a resonant peak for the q ^2 
mode at 0.663 per rev with an associated damping factor 0.091. 

The actual values in the simulation for the q ^2 modal frequency and 
damping were 0.666 per rev and 0.043, respectively. 

The maximum likelihood method was applied to the same data. 

This identification method resulted in an estimated q^2 
quency of 0.671 per rev and damping of 0.036. These estimates 
compare favorably to the true values used in the simulation. 
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Figure 5.1 shows the time histories for the collective pitch 
random input and the wing chordwise acceleration measurement from 
the nine degree-of- freedom simulated data. Superimposed on the q ^2 
measurement is the estimated measurement time history from the 
simple second order model. 

Since no data windows were used for the time series analysis 
method, it should be concluded that the apparent discrepancy of 
the result (relative to the simulated and maximum likelihood esti- 
mates) could be reduced significantly. 

In order to evaluate the maximum likelihood method further, 
the second set of data (using sum of sines inputs) were processed. 
The results were an estimated niodal frequency of 0.663 per rev 
and damping of 0.033. These results and those of the spectral 
analysis are summarized in Table S.l. 


Table 5.1 

Frequency and Damping Estimates of q ^2 Mode 



03 


Values used 1n simulation 

0.666 

0.043 

Values from spectral analysis of 
random input case 

0.663 

0.091 

ML values identified from random 
input case 

0.671 

0.036 

ML values identified from sum of 5 
sines input case 

0.663 

0.033 


5.4 CONCLUSIONS 

This brief example has demonstrated: 

(1) the possibility and usefulness of generating simulated 
wind tunnel data as a data base foi which the true 
values of modal frequency and damping are known. 
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(2) a procedure for assessing the accuracy of methods (e.g., 
spectral analysis or parameter identification) for esti- 
mating modal frequencies and damping factors, and 

(3) the ability to assess which types of control inputs will 
result in better estimates (random and sinusoidal inputs 
were used above) . 

Furthermore, the introduction of aerodynamic gust effects to the 
simulation, which was demonstrated in Chapter II, adds the possibil- 
ity of studying their effects on the results of various data reduc- 
tion techniques. 
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VI. RECOMMENDATIONS 


This study has achieved the following objectives. First , 
several advanced analytical tools have been successfully demonstrat- 
ed on the XV- IS mathematical model, including hijh order transfer 
function evaluation and maximum likelihood parameter identification 
techniques. Second , the application of this analysis has produced - 
a preliminary guideline for selecting test inputs and instrumenta- 
tion. Third , specific problems have been isolated which may limit 
the information which can be extracted from the tunnel tests. 

It is recommended that several aspects of this preliminary 
study be expanded. These include the following: 

(Ij Incorporation of existing computer programs to determine 
the sensitivity of transfer functions Ce-g«» frequency 
response) to particular parameters such as tunnel velocity, 
wing stiffness, collective pitch, or measurement error. 

This would produce a valuable guide to estimating the 
most significant error sources in determination of fre- 
quency or damping. 

(2) Investigation of the effects of nonlinearities in wing 
structural parameters, the control system, or aerodynam- 
ics. These nonlinearities would include backlash, hys- 
teresis, or deadband. 

(3) Determination of effects of fuselage modes on wing 
response. The present support admittedly limits the 
effect of fuselage coupling, but these couplings may be 
important in analyses of flight data. 
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One further reconmendation is for the development of an ad- 
vanced technique for processing data to determine, on-line, fre- 
quency and damping. This recommendation is based on the results 
of Chapters III to IV of this report, and evaluation of existing 
methods of spectral analysis. The basic requirement is to develop 
the capability to determine an estimate of the required analytical 
model structure and parameters which best approximate a multivari- 
able modal response. A fundamental result of this report, for 
example, is that multivariable response poses a difficult problem 
in establishing modal contributions from a particular transfer 
function. This problem will be further magnified by the require- 
ment to determine frequency and damping from multivariable response 
data. 


The following procedure is, therefore, suggested to implement 
this recommendation: 

Cl) Develop algorithms which provide on-line estimates of 

model structure Ce*g*» transfer function) and parameters 
of that model. 

(2) Evaluate this algorithm on the simulation discussed in 
Chapter II, and the transfer function evaluation of 
Appendix B. 

C3) Implement the algorithm on-line at Ames Research Center 
and evaluate it on helicopter model data. 

(4) Perform any modifications or further extensions of the 
algorithm. 

It is anticipated that successful development of such a tech- 
nique, combined with the results of this report, will produce a 
state-of-the-art technology for improving the results of advanced 
rotorcraft tunnel testing. 
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APPENDIX A 

SUPPORT EQUATIONS OF MOTION 


The addition of the support degrees of freedom to the equations 
of motion is accomplished in a manner very similar to the inclusion 
of the wing degrees of freedom as described in Ref. 3. The equa- 
tions of motion of the support and wing will be found as functions - 
of the forces and moments acting on the hub due to the rotor. 

These can then be combined with existing equations (in Ref. 3) of 
motion for the rotor, which require the motion of the hub, to com- 
plete the equations. 

The equations for the wing and support degrees of freedom may 
be found by use of Lagrange's equation 



where 

L - T-V 

T ■ Kinetic energy of the system 

V ■ Potential energy of the system 

■ Generalized force for i^ generalized variable 
generalized coordinate 

Computing the kinetic energy of the system requires the iner- 
tial linear and rotational velocities of the pylon center of mass 
and the fuselage center of mass, which are computed below. 

Figure A.l shows the geometry and coordinates of the wing anc 
fuselage system. The z^ axis is the rotational axis of the 'i'gyp 
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Figure A.l Support and Cantilever Wing Geometry 
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degree of freedom. The fusealage centsr of gravity is a distance 
Ixf forward of this axis. The cantilever root restraint of th." 
wing elastic axis is at point A, which is forward and 

laterally (posiuive right) from the axis. The wing has sweep 
angle 6^3 (positive aft), dihedral 6 ^^ (positive up), and in- 
cidence angle 6vf (positive leading edge up). The wing semi- 
span is /T;^* wing tip the '^W]^ » ^W2» ^W3 angles are 

reversed so that the rotor hub axis, z^, is parallel to x^ . 

The rotor forces and moments acting at the hub are T, H, Y, Q, 

M , M , as shown. 

X y 

The system degrees-of-freedom are. 


X ■ translation in x_ direction 
s o 

y^ = translation in y^ direction 

4; « rotation about z axis 

sup u 

■ lowest mode wing vertical bending (positive up) 

qw2 “ lowest mode wing chbrdwise bending (positive up) 

p • lo.vcot mode wing torsion (positive leading edge 

up) 


The inertial reference frame is taken to be colinear with the 
frame centered at point 0 as shown in Figure A.l and fixed in 
inertial space. 


Given the position of 
noted X, its position in 


a point expressed in the A frame, 
the inertial frame is: 


de- 
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/ 
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c-i»3 
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c-^s 
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(A. 2 ) 
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Similarly, given an angular position vector expressed in the A 

frame, denoted a, it is expressed in the inertial frame as 
A 


a 

I 



(A. 3) 


Using these two equations and equations in Refs. 1 and 2, which 
express the position of the hub, pylon c.g., and an arbitrary ele- 
ment of mass in the wing with respect to the A frame, it is 
straightforward to express these positions in the inertial frame. 


The task of computing the kinetic energy will be much simpler 

• T * 

if XX is first calculated for an arbitrary x. Later, the pylon 
I I A 

and wing positions, x_ and x , will be substituted for x. 

aP a'' a 

T 

Let X “ [x, y, z], then from Eq. (A.l) 

A 


X 

I 




CA.4) 


X 

I 


" Xs+zcil^g-ysij/s-ig ( (2+«>Xw) si's* (y+^y„) ct|>g) 

ys+2si|»g+ya*.s*’t's CC2+ixw^<^'»'s’ 

- -X 


(A, 5) 
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The position of the pylon center of mass in the A frame is: 



; 

■i 


f 




K 


! 

[ 


yVwi+qwi Cy-hpn63)+qw2 (‘y<S2*hp^'Si) 

* P„[V(l'EA-''p)**EA^2l 


yT„*1wi [‘hpTS 2 -yS ^ 1 [hpn -/« 3 I *P„[-hp« i»hp« 3 « 2 1 
-yT„« 3-<lwiy« 2 •‘'W2>'*Pw t‘ ( t'’EA''’p) ^ 2**EA> ' *''p 

(A. 6) ■ 


where the shorthand notation of Ref. 3 has been adopted. Similar- 
ly » 


a 

I 


P 


■-nq^ + ri62^2 " *^3^ 

-'"sup * '’'* 2 'll * '"*2 



(A. 7) 


To compute the kinetic energy of the wing, the velocity of a 
differential mass element will be found, and the result integrated 
over the wing. Let this differential element be a distance z 
ahead of the elastic axis and a distance r from the root measured 
along the elastic axis. Then the position of this element in the 
A frame is: 


f " 


! ■ 

I . 


: . 


I 



r6wi+25w2+qwiCn(r)-zn'(r)63)+qw2(-n(r)62'^zn'(r)5]_) 

+PwC5Cr)-C^(yB^)r63) 


r+q^^Czn" Cr)52*qC’^)^l)^‘^W2(2’l"C^)‘^Cr)63)-p^5(r)6j^+z5 3 


-r6w3+z*qwinCr)6 2-qw2'lC’^)'Pw^ 


(A. 8 ) 
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Finally, the position of the fuselage in the inertial frame is 



(A. 9) 


(A. 10) 


All the parts are now available to compute the kinetic energy. 
The positions in the A frame must be expressed in the inertial frame 
through Eqs. (A. 2) and CA.3). Then the positions must be differen- 
tiated with respect to time, and the resulting velocities substi- 
tuted into the kinetic energy of the system, which is 


T = 1 m 4. ^ TH *1 Ty .1 Tt 

^ I P^P^P I * J Ip“p * I 



•T* 

WWW 


dz dr 


The potential energy is simply 


I * >'q2'l2 * * Kxs*, * hsA * 


+ K(l) i|) 
^sup^sup 


(A.ll) 


(A. 12) 


The generalized forces include aerodynamic forces on the wing and 
fuselage and structural damping terms. The aerodynamic, forces 
for the wing were derived in Ref. 3. The aero forces on the 
fuselage will be neglected because they are small by comparison 
to the wing and rotor forces. 
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By substituting Eqs. (A. 11) and (A. 12) into Eq. (A. 1) and 
performing the indicated differentiation and by neglecting small 
terms, the following vector equation results: 


“ ^""^f (A. 13) 

where 

x„ - [q«i q «2 p« 1-sup *3 
5£ ■ wing flaperon deflection 

longitudinal gust 
lateral gust 
vertical gust 

2Cjj 2Cy Cq 

3x4 matrix from p ... of Ref. 3. 

3x4 matrix 



[ bgl 3x3 matrix from p ... of Ref. 3. 

7” J 3x3 matrix 

and the matrices a 2 » aj^, a^, and a are given on the following 
pages. 


I : 
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Si. 


, is total yaw moment of the aircraft 

exclusive of the rotors. 


m. 


w 


Si. 


m is mass of the wing 
w 




VTw_, m£ is mass of the fusealage 


Si. 


*’Yw " ^'yw^^Tvf* ^Yw lateral component of distance 

from wing root elastic axis to c.g., 


'Xvr * ^Xv^^Tyf* ®-Xw longitudinal component of distance 

from wing root elastic axis to aircraft. 


c.g. , 


and all other parameters are as defined in Ref. ‘3. 

The total equations of motion are found by combining the 
rotor equations and Eq. (A. 133 as described in Ref. 3. 

The numerical results of Chapter II were generated with the 
following simplifications: 

Ca) ® torsion of supports) 

(b) ■ 0 (wing cantilever attachement at c.g.) 
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APPENDIX B 

TRANSFER FUNCTION ANALYSIS 

B.l REVIEW OF MATHEMATICAL FORMULATION 

As explained in Chapter II, the mathematical model for the 
XV-IS, as developed by Johnson [2-5], expressed in vector notation,, 
is 


X - Fx + Gu + rv 

■ Hx ♦ Du (measurements excited by the control) 

» Hx ♦ D'y (measurements excited by random gust) 

There is an F, G, V , D, and D' associated with the symmetric 
motion of the rotor/cantilever wing, and another F, G, T, D, and 
D' associated with asymmetric motion of the rotor/cantilever wing. 
The symmetric and asymmetric models are combined (as explained in 
Section 2.3) for a complete description of the dynamics at the left 
and right side of the aircraft. 


(B.l) 


(B.2) 

were calculated using the Leverricr method (9,10] to evaluate the 
adjoint. The transfer function was the primary analytical tool 


The frequency responses 


and 


I(s) 
uTsT ’ 

rcs) 

y(s) 


H(sl-F)’^G + 


H(sl-F)‘^r 


+ D' 


(referred to as signal 
transfer functions) 


(referred to as noise 
transfer functions) 



1 


I 


used to evaluate the dynamics of thj XV- IS and» hence, evaluate the 
instrumentation requirements and test guidelines for the wind tunnel 
testing of the XV-15 tilt rotor aircraft. 

B.2 DERIVATION OF FORMULAE 

Transfer functions for measurement of the states or linear 
coordinates of the states were used in several ways : 

(1) The transfer function indicates if a particular measure- 
ment excited by a particular control shows a resonant 
peak to identify a desired mode. 

(2) To study the stability of the XV-15, estimates of the 
damping (c) of the least damped modes are desired. For- 
mulae used for calculating c from transfer function data 
are discussed in Section B.2.1. 

(3) To quantify the value of different transfer functions 
for identifying the frequency and damping of a desired 
mode, and to clarify which modes contribute substantially 
to a resonant peak (when more than one mode are close 

to the peak frequency) , the modal power and modal power 
ratio were derived. These formulae are presented in 
Section B.2. 2. 

(4) To analyze the effect of various gust levels in the wind 
tunnel on the ability to identify the frequency and damp- 
ing of a desired mode, the noise-to-signal ratio was 
derived as explained in Section B.2. 3. The signal power 
and noise power were calculated from the signal and noise 
transfer functions over the same frequency range. 
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B.2.1 Damping Calculation [6] 

To evaluate the damping of a particular mode, the transfer 
function is approximated by a general second>order system in a 
region about the resonant peak for that mode. The second order 
system is 


mx + cx kx ■ f 


CB.3) 


where 


0)^ ■ i/Tc7m is the natural frequency 


5 - c/2Wj^m is the damping factor 


f ■ control input. 


This equation may be rewritten as 


mx ♦ cx + mw^x ■ f 
n 


(B.4) 


The transfer function from f to the system acceleration, a ■ 'f, is 


Iff} - T(s) 


J 

ms + cs + mw„ 
n 


(B.5) 


The Bode magnitude plot requires the magnitude of this func- 
tion evaluated on the ju axis of the s-plane. 


T(ju)) 


-« 


(B.6) 


m(a)^-w ) ♦ jcu 


1 
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TCjca) 




CB.7) 


Eqs. (B.6) and (6.7) may be combined to yield 


T(jw) - - [m(«^-w^) - jcwl 

/ti • 


(B.8) 


This form is preferred over Eq. (B.6) because the complex denomina< 
tor is eliminated. 

Taking the imaginary parts of both sides of Eq. (B.8), inte- 
grating both sides, and solving for the average value of c gives 




[' * 3 ^ 


(B.9) 


where the interval (oi) 2 »u 2 ) is chosen such that <d^ is at its midpoint 
and for which the approximation of the frequency response curve 
by the second-order system is valid (e.g., ■ 0.8 and U 2 ■ 

1.2 0 )^). The frequency at the peak, Up, is used as the estimate of 
which is reasonable for small values of ; since 


«i»_ ■ w_ /1-2C* 

P n 

for a second order system. 


(B.IO) 
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It is now necessary to compute an average value for m in order 
to compute the damping factor from 


2b)_m 

II 


(B.ll) 


To this end, Eq. (B.S) is written as 


T(s) - ^ 


s^/m 


5 ♦ 2c«_s ♦ tt) 

n n 


(B.12) 


This leads to an equation analogous to Eq. (B.7); that is, 


TCju) I ^ ■ "T 


U) 


+ C2C(i)j^w) ^ 


(B.13) 


Substituting 2cu^m for c in Eq. (B.S) gives 


T(j(u) ■ - fmCw^-w^) - j25u^mwl 

(l) ^ ^ 


(B.14) 


Substituing Eq. (B.13) into (B.14) for |T(jw)r and then taking the 
imaginary parts of both sides of Eq. (B.14), integrating both sides, 
and solving for the average value of ra gives 


m 


(I). 

/: 


I»(T) 


(1) 


du 


(B.15) 


where 


u). 


■ nc) ■ / — 2 TT 


du) 


(B.16) 
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This integral may be computed analytically as follows: 




U)_ -« 

1 . n 


tan-l “4 


2 /i^C * 


Finally^ substituting Eqs. (B.d) and (B.15) into Eq. fB.ll) results 
in an implicit equation for the damping. 


(B.17) 


or 

C - (B.18) 

where P is defined appropriately from Eq. (B.17). P is a function’ 
of the transfer function between the integration limits of and 
b> 2 * and is not a function of Therefore, c may be found by 
iterative use of Eq. (B.18). It was found that, given a good 
initial guess for convergence was achieved in three of four 

iterations. 

A good initial value for ; is 

'o ■ TO ■ ro V (?) ■ 

which was found from Eq. (B.18) by expanding the analytic e.xpression 
for 1(c) in a power series and neglecuxng second and higher order 
terms in c* 
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The transfer function T(ju») was approximated by a second 
order system in the region of the resonance. A second order system 
has a phase angle of -90® at the peak. Thus, the effect of other 
modes on the resonance can be lessened by shifting the phase of 
the measured transfer function before applying Eq. CB.17). This 
phase shift is performed by the following equation: 


T * T 
new 


* 

-•5 P 


CB.20) 


where Tp is the measured response at the resonant peak. 


Data Processing 

The calculated value of 5 is a function of bandwidth since 
other modes interfere as the bandwidth increases and the accuracy 
of the integrals (averages) used in the calculations may deteriorate 
as bandwidth decreases, particularly if the number of data points 
in the bandwidth is small. To investigate this dependency, the 
following procedure was used. 


Using the transfer function as an example, T(ju)) was 

computed for 100 values of oj between 0.3S and 0.45 per rev for the 
symmetric case (qj^^ is 0.398 per rev), and 100 values of w between 
0.49 and 0.62 for the asymmetric case (q^^^ is 0.562 per rev). Using 
these data, 5 was calculated for various bandwidths containing each 
peak. The results are summarized in Table B.l. 


The calculated values of damping presented in the body of this 
report were based typically on 10 to 15 data points in the frequency 
range of integration (i.e., in the bandwidth purpose 

of obtaining 100 points in this small range was to reduce the ntimer- 
ical errors in the calculation to a negligible level. Therefore, 
the error in the calculated values of ? are due to the inteiference 
of other modes and to the approximations in the theoretical method 
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Table B.l 

Calculated Damping as a Function of Frequency Bandwidth 

Used for Integration 


BANDWIDTH 


5 FOR MODE FROM 
(iwi/Sf TRANSFER FUNCTION 


•-I 

“2 

SYMMETRIC MOTION 

(D„ * 0.4002 
P 

ASYMMETRIC MOTION 
(Dp = 0.5616 

0.99(d„ 

p 

- l.OKDp 

0.04098 

0.04200 

0.98(d„ 

P 

- 1.02(d^ 
P 

0.04331 

0.04619 

0.97(d„ 

p 

O. 

o 

« 

t 

0.044B6 

0.04977 

0.96(d„ 

P 

- 1.04(Dp 

0.04567 

0.05194 

0.95(Dp 

- 1.05(d 

0.04584* 

0.05302 

0.94(Dp 

CL 

O 

* 

1 

0.04578 

0.05348 

0.93(Dp 

- 1.07(Dp 

0.04503 

0.05361* 

0.92tDp 

- 1.08(d„ 
p 

0.04537 

0.05359 

0.91(Dp 

- 1.09(Dp 

0.04511 

0.05349 

O.OODp 

- 1.10(Dp 

0.04483 

0.05336 

True Value of c 

0.04800 

0.05540 

% Error (Based on 
Best Value of ;) 

5.0% 

3.2% 


Maximum value of ; 


of the calculation. Further study of data processing techniques 
should include investigation of other methods, such as that suggested 
by Kenney and Pancu [11]. 

It was found that with 10 to IS data points in the region 
0.8 (Dp to 1.2 (Dp that 0.9 (Dp to 1.1 (Dp gave the best value for damp- 
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ing and with 100 data points in the 0.8 Up to 1.2 Wp band, that 
0.95 Wp to 1.05 Wp gave the best value for damping. Given a small 
frequency increment between data points, the calculated damping does 
not vary significantly in the range 0.9 Up - 1.1 Wp to 0.95 Wp - 
1.05 Up . 

This analysis was performed on the (0,398 per rev) and q^^^ 
(0.562 per rev) modes. The trends shown in Table B.l may change 
somewhat for high frequency modes (such as the mode at 2.43 per 
rev), where the same frequency ratio about Up gives a bandwidth 
several times greater. 

B.2.2 Modal Power and Modal Power Ratio 

To best identify particular modes from the frequency responses 
of the system, it was desired to: 

• Quantify the quality of different transfer functions to 
identify a particular mode 

• Quantify the extent to which a mode is obscured when more 
than one are contributing to a resonant peak on the fre- 
quency response. 

The quantification of a particular mode’s contribution to the 
transfer function can be accomplished in the following manner. 
Consider the response of a transfer function to sinusoidal excita- 
tion and express this response in terms of the response to an im- 
pulse function. In other words, express the steady state response 
to a sinusoidal forcing function in terms of transient response 
quantities, namely the residues at the poles of the particular 
transfer function. 
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The response to these two inputs will be expanded in a partial 
fraction expansion. Let the subscripts t denote transient, associ- 
ated withan impulse input; s denote sinusoidal, associated with a 
sinusoidal input; and ss denote steady state, associated with a 
sinusoidal input after the transients have died out. 

For an nth order system, 

y^Cs) » T(s)Z{5(t)} » T(s) 5 (B.21) _ 

n k. 

YtCs) - Z txttT 
^ i»l ^i^ 


where N(s) and D(s) are the numerator and denominator polynomials, 
respectively, and the X^'s are the poles of TCs) and the k^'s are 
the residues at those poles. 


Forcing the system with a complex sinusoid of unit magnitude 


gives 


y^Cs) 


T(s) 


Ts'-fS;:) ■ 


(B.23) 






. ss 


(B.24) 


From Eqs. (B.21) and (B.23), it is seen that 
y^(s) - (s-ju)Q)yg(s) ■ 0 


(B.25) 
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The next step is to substitute Eq. (B.22) for y and Eq. (B.24) 
for in Eq. CB.2S) and to express the result as the ratio of two 
polynomials, whose numerator polynomial must be *0 identically 
for all values of s. Hence, the coefficient of each power of s 
must » 0. The following relationships result. 


n 

k__ - Z k: 
i-1 


.r- 




Therefore, 


n k. 


w * • Z g r 

ss 


The steady state response is then 


^ssf"^ " ^ss® 


jtDot , gjwot 2 *^i 


i-1 


(B.263 


(B.27) 


(B.28) 


(B.29) 


A complex sinusoid was used because it simplifies the algebra. 
The steady state response to a purely sinusoidal excitation, 
u • sin w-t, can be obtained by merely taking Im[y (s)]. 

The magnitude of the steady state response is 

|T(jo,^)l - kg3 (B.30) 

where TCs) is the system transfer function. The power in an 
interval (u)2»Wj^) about a resonant peak is defined as 
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P M iTCjw)!^ do) 

“2 

- j f(jw) • f(j(u) du) (B.31) 

“1 

where f(ju) is a vector in the complex plane and CO indicates the 
vector inner (.dot) product. 

From Eqs. (B.29) and (B.30), 


where 



(B.32) 


(B.33) 


This is equivalent to saying that the total system response at 
frequency u is the sum of the responses of the individual modes 
(i.e., X^’s). 


For each oscillatory mode, there exists a conjugate complex 
pair of poles, which will be denoted X^ and » X^, whose associ- 
ated residues are k. and k. ■ Ic . . Let 

X XX 





Tj^-jw 


(B.34) 
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Then the contribution of this oscillatory mode to the total power 
is 


‘’i 





T(jw) du 


CB.3S) 


and ratio of this fractional power to the total power (Eq. (B.31)) - 
is called the modal power ratio for the ith oscillatory mode; 

Pi 

MPR - (B.36) 

B.2.3 Noise-to-Signal Ratio 

The transfer functions given in Section B.l 

•*- H^Cjw) (signal transfer function) 

-► Hj^(ju)) (noise transfer function) 

give the gain and phase of the system when excited at a discrete 
frequency, and since the model is linear, any number of frequency 
responses can be added to give the response to the sum of the 
inputs. Theoretically, white noise would give the entire transfer 
function since it includes all frequencies. 

Thus, to establish the ratio of how a particular measurement 
responds to random gusts in the tunnel compared to how it responds 
to commanded controls. The signal squared, or power of the two 
transfer functions, can be compared in the frequency range of 
interest. 


lOS 



i 


f 


r^2 

I |u|^ |H. 


do) 


S/N 


r 


q dw + r((ij- • w, ) 


Wi 


where q is the gust covariance and r is the measurement covariance 
(considered as white noise in the frequency range of interest) 

q • 2 T x(t) 

^ corr ■' 

where 

Tcorr correlation time of the noise ■ L/u^ 

(correlation distance divided by the wind 
velocity) 

x(t) is the autocorrelation of the noise 

X(0) - ''Ims 


is a composite power spectrum obtained in the follow- 
ing way: 


■ <j>l 

out -^in 

iHjI 

" ^ 2 . 

out in 

iHjI 

out •’in 

IH 3 I 

correspond 

to t 


vertical spectrum, respectively, of the Dryden model, and H 2 
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Hj are the gust transfer functions from these three sources to the 
desired measurement. Assuming the gust from the three directions 
is uncorrelated means that the power from these three signals can 
be added in an RMS sense; therefore, these sums can be integrated 


<t> 


n 


3 

Z 

i-1 


<^i 

out 


S/N 


u 


rms 



|H|^ dw 



(j»^ do) + r(u)2-Wj^) 


(B.37) 


The reciprocal of this relationship, the noise-to-signal ratio 
(NSR) , is more convenient to work with since the gust effects can 
be studied alone without the measurement noise. 


N/S 


rms s 


r(u)2 

T" 


■Wj) 


“rms^s 


(B.38) 


where and are the power under the noise and signal power 
spectra. The noise>to-signal ratio for the measurement can be 
added in if it's characteristics are known. It should be pointed 
out that this is the noise-to-signal ratio for one data record. 
Normally, k records are recorded and averaged in the frequency 
domain, where the signal adds and the random noise tends to cancel, 
greatly improving the signal- to-noise ratio. Bendat and Piersol [12] 
suggest using 10 records or more. The principal restriction is the 
total sample time. 
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The point to be emphasized is that the signal-to-noise ratio 
in Eq. (10) is for one record, and thus all the noise -to-signal 
(NSR) plots in Chapters III and IV are for one data record. 
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APPENDIX C 

TRANSFER FUNCTION FREQUENCY RESPONSES 

This appendix presents the Bode plots for the rotor/cantilever 
wing and coupled wing models. These plots are outputed from the 
tr{ nsfer function analysis program. Other data outputed are the 
poles, zeroes, and residues of these transfer functions. These 
latter outputs are not presented here. 

Each Bode plot is given with a reference scale which indicates, 
according to the poles of the transfer function, which degree of 
freedom is predominant at a particular frequency. 


109 































































































Figure C.29 “ Wing Vertical Acceleration to Collective 

(Coupled Wing Response • Excite One Wing, Measure 
Same Wing) 
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Asymmetric: 


Symmetric: 
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Figure C.30 Chordwise Acceleration to Collective 

(Coupled Wing Response * Excite One Wing, Measure 
Same Wing) 
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Figure C.31 * ^ing Vertical Acceleration to Flaperon 

(Coupled Wing Response * Excite One Wing, Measure 
Same Wing) 
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